The disordered flat phase of a crystal surface- critical and dynamic properties 
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We analyze a restricted SOS model on a square lattice with nearest and next -nearest neighbor 
interactions, using Monte Carlo techniques. In particular, the critical exponents at the preroughen- 
ing transition between the flat and disordered flat (DOF) phases are conflrmed to be non-universal. 
Moreover, in the DOF phase, the equilibration of various profiles imprinted on the crystal surface is 
simulated, applying evaporation kinetics and surface diffusion. Similarities to and deviations from 
related findings in the fiat and rough phases are discussed. 
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A new type of crystal surface, the disordered flat 
(DOF) phase, has been introduced by den Nijs and Rom- 
melse almost twenty years ago.^'^''^ It is characterized by 
finite height fluctuations, although it comprises a dis- 
ordered array of steps, thereby combining features of 
the, conventional, flat and rough phases. Experimen- 
tal evidence for a DOF phase has been reported, for in- 
stance, for films of rare gas atoms on graphite and for 
GaAs(001)ii^ 

A minimal model to describe the DOF phase has been 
proposed by Rommelse and den Nijs.^ It is a restricted 
SOS, or RSOS, model on a square lattice with interac- 
tions connecting pairs of nearest and next-nearest 
(j, k) lattice sites, with the Hamiltonian 



0.5 







1 


reconstr.\ 
rough \ 

- -V 


rough 
•' 


*' 


V 

y 

reconstructed 
flat 

.1.1. 


DOF .4** 
T 

. •/ . 1 


flat 



(1) 



where hi is the integer- valued height at site i. The ab- 
solute value of the height difference is restricted to be 
at most one for neighboring sites, and, thence, at most 
two for next-nearest neighbors. 6 is the Kronecker sym- 
bol. The interaction between next-nearest neighbors, L, 
is assumed to be positive. 

The phase diagram of the Hamiltonian (1) in the 
{kBT/L,K/L)-pla.ne is depicted in Fig.l. It displays 
five distinct phases: The usual flat and rough phases, 
the reconstructed flat and rough phases arising from a 
repulsive, negative nearest-neighbor interaction, K < 0, 
as well as, finally, the DOF phase. The most intrigu- 
ing aspects of the phase diagram are the disordered fiat 
phase and the transition between the fiat and the DOF 
phases, the preroughening transition^ 

Similar phase diagrams have been obtained for re- 
lated SOS as well as other models with short-range 
interactiongSiiiS, including variants to describe fcc(llO) 
surfaces. In the following, we shall first discuss critical 
properties at the preroughening transition and then dy- 
namic properties of the DOF phase, analysing the min- 
imal model (1). The findings are expected to be largely 
independent of details of the model. 

To determine critical properties at the preroughening 
transition, we used Monte Carlo techniques, in particu- 
lar, apart from the standard Metropolis approach, the 
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FIG. 1: Phase diagram of the RSOS model, eq. (1), as 
obtained from a transcription of the diagram determined by 
Rommelse and den Nijs-, with the circles referring to their 
results using the transfer matrix method. The boundary lines 
of the DOF phase at low temperatures follow from the present 
Monte Carlo simulations. 



Wang-Landau and the transition matrix algorithmsiSiili 
Square systems of size Af^, with M up to 64, have been 
simulatediii Taking into account finite-size effects, we 
estimated standard critical exponents, namely those of 
the order parameter (the parity?), /3, and its susceptibil- 
ity, 7, as well as of the correlation length, i/, at K/L~ 
1/8, 1/4, and 1/2, see Fig. 1. 

Note that the preroughening transition temperature 
Tp, at fixed ratio r = K/L, may be obtained quite ac- 
curately by fitting the height fiuctuations to the ansatz 
Ceff ln(M). Increasing M, Ceff tends to zero in the fiat 
and DOF phases, while it seems to approach at Tp a non- 
vanishing value, with only minor finite-size dependences. 
Thence, the height fluctuations seem to diverge directly 
at the preroughening. The prefactor Ce//(M — > oo, Tp) 
is estimated to be smaller than the characteristic value 
I/tt^ of a transition of Kosterlitz-Thouless type (which 
occurs, e.g., when going directly from the flat to the 
rough phase, see Fig. 1).^^. It is found to increase with 
K/L, reaching, possibly, I/tt^ at the triple point between 
the flat, DOF, and rough phases. 
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The critical exponent (3 is found to tend to increase 
with the ratio of the coupUngs r — K/L from 1.08 ±0.09 
at r=l/8 to 1.16±0.02 at r= 1/4, and then to 1.75±0, 03 
at r =1/2. The corresponding estimates for 7 are 
2.75 ± 0.23, 2.80 ± 0.03, and 3.50 ± 0.05. For v we 
obtain 2.40 ± 0.20, 2.52 ± 0.03, and 3.43 ± 0.05. The 
rather unusual values of the critical exponents at the 
preroughening transition are in accordance with previous 
estimatesSi^. Our findings confirm the non-universality 
of the preroughening transition^, with the critical expo- 
nents depending continuously on K/L. Applying scaling 
laws, the critical exponent of the specific heat, a, turns 
out to be strongly negative. Actually, from our simula- 
tional data, the specific heat appears to vary smoothly 
with temperature, preventing a direct estimate of a. 

Dynamic properties of the DOF phase have not been 
studied in much detail before, despite a wealth of theoret- 
ical and experimental analyses of morphological changes 
for flat and rough crystal surfaces. In fact, we are aware 
of only two studies on the equilibration of initially flat 
surfaces in the DOF phase, driven by evaporation kinet- 
ics, using either Monte Carlo techniques (studying the 
Hamiltonian (1) with K = Q) or the renormalization 
group methodii2ii2. 

In these two studies, the time evolution of the height 
fluctuations < 5h? > has been analyzed (iii^ For systems 
of A'P sites, one defines < Sh'^ >=< Y,{h, - h*f/M'^ >, 

i 

where h* denotes the average height. The height fluc- 
tuations have been observed and arguediSii^ to increase 
initially logarithmically in time, both in the rough and 
DOF phases of RSOS models. In the DOF phase, then 
the fluctuations decrease, and finally they approach the 
equilibrium value. In contrast, in the rough phase, the 
fluctuations grow monotonically before they eventually 
saturate, for finite M . 

We confirmed these findings, for K = and, in the 
DOF phase, ksT/L — 0.7 (2.5 in the rough phase), with 
M ranging from 16 to 256, using full periodic bound- 
ary conditions (previously helical boundary conditions 
had been employed)^ The initial height configuration 
is hi{t = 0)=0. To mimic evaporation kinetics, we ap- 
ply the standard Metropolis Monte Carlo algorithm, at- 
tempting in each move to change, at a randomly chosen 
site i, the height by ±1. In addition to the previous work, 
we identified, at very early times, a diffusive behavior, 
with < 5h? >oc t, due to uncorrelated height changes at 
randomly chosen sites. Of course, it seems doubtful that 
the logarithmic behavior of < 5h^ > (t) will prevail to 
arbitrarily long times in the DOF phase when studying 
larger and larger systems, as it may be the case in the 
rough phase. 

The time evolution of the surface may be visualized 
by monitoring Monte Carlo configurations. In the DOF 
phase, at early times, several clusters of different lo- 
cal average height, 1/2 or -1/2, are formed, with the 
height fluctuations growing logarithmically. In the clus- 
ters, the sites are, in a disordered fashion, predominantly 
at heights and 1 or and -1. The subsequent decrease 




FIG. 2: Typical equilibrium Monte Carlo configuration in 
the DOF phase, with K = and ksT/L = 0.1, comprising 
sites at height -1 (black), (dark grey), and 1 (light grey). 
The surface has 24^ sites. 

in < Sh'^ > is due to the fact that one type of cluster, 1/2 
or -1/2, prevails in the equihbrated DOF phase, see Fig. 
2. Of course, with full periodic boundary conditions, the 
overall average height may eventually wander, tending 
to stick at neighboring half-integer values. This effect, 
with the entire surface undergoing a random walk, does 
not affect the height fluctuations. 

To analyze more completely the dynamics in the DOF 
phase, we also studied the equilibration of various other 
one-dimensional proflles imprinted on the surface, in 
particular, steps, periodic gratings, and wires^ There, 
the interplay of the orientation dependent surface ten- 
sion and mobilitjii^ii^ leads to interesting morphological 
changes in the rough and flat phases, for evaporation 
kinetics as well as for surface diffusion, see pertinent re- 
views on theory and experimentssiSiiLiSiiS 

In the DOF phase, a step of monoatomic height may 
be introduced by fixing the average heights of the sites 
at two opposite boundaries to be, e.g., at 1/2 or 3/2, by 
setting the heights at these boundary sites alternately 
and 1, or 1 and 2, respectively. The other two boundaries 
are connected by periodic boundary conditions. Square 
systems with sites are considered. Averaging over 
sites parallel to the fixed boundaries as well as over many, 
N, realizations of the equilibration process, one obtains 
a one-dimensional height or step profile H{x,t), with 
X running from, say, -M/2 + 1/2 to M/2 - 1/2. At 
time t = 0, the height configuration is assumed to look 
like a checkerboard comprising heights and 1 in the 
half of the surface next to the 1/2-boundary, and heights 
1 and 2 in the other half. Then, H{x,t = 0) = 1/2 
for a; < 0, and H{x,t = 0)= 3/2 for x > 0, with a 
sharp step in between. As time proceeds, the step profile 
will start to broaden, see Fig. 3. In the flat and rough 
phases, the step width is knowni^ to grow asymptotically 
in time, t and in the thermodynamic limit, M — > 00, 
with a power law oc t^. The characteristic exponent b 
depends on the transport mechanism and the type of 
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FIG. 3: Evolution of step profiles, H{x, i), in the DOF phase, 
K/L = 0, ksT/L = 0.1, at various times (measured in Monte 
Carlo steps per site, MCS), applying evaporation kinetics. 
Systems of 128^ sites have been simulated, averaging over 
A'^ = 11200 realizations. For clarity, only each third profile 
point, on both sides of the step, is shown. 



phase. For instance, in the case of evaporation kinetics, 
6 is 1/2 in the rough and 1/4 in the flat phase, in the case 
of surface diffusion, b is 1/4 in the rough and 1/6 in the 
flat phase iSSiSiiSSiSiSi. b may be determined from Monte 
Carlo data by monitoring the position of the step profile 
at fixed height iS^ In that way, an effective exponent beff 
may be estimated, approaching the true value b in the 
limits M, t — > oo. 

Typical simulated step profiles, applying evaporation 
kinetics, in the DOF phase are depicted in Fig. 3. 
Analysing such profiles, wc obtain an effective exponent, 
for M = 128, approaching, from above, closely 1/2 in 
the rough phase. In the DOF phase (M = 128, K = 
0, fc^r/L = 0.1), the effective exponent also decreases 
monotonically with time, the data being well compatible 
with an asymptotic value b = 1/4, as in the flat phase. 
Surface diffusion may be mimicked in the simulations by 
attempting to, say, decrease the height at a randomly 
chosen site by one and to increase simultaneously the 
height at a neighboring site by one. For that dynamics, 
the effective b is found to approximate nicely 1/4 in the 
rough phase, M = 128, and to approach 1/6 in the DOF 
phase, M = 256, K = 0, ksT/L = 0.1. 

An interesting phenomenon occurs when periodic, one- 
dimcnsional gratings arc imprinted on the surface. In the 
flat phase, the healing proceeds layer by layer, driven by 
an annihilation of the meandering steps bordering the 
top and bottom terraces of the gratingi ^yi^^i^^i^Ti^?!^? . 
To study the equilibration of gratings, especially in the 
DOF phase, we consider an initial profile of the form 
H{x,t = 0) = aint{Asm{2TTx/M)), i.e., a discretized 
sine function with amplitude A. Full periodic boundary 
conditions are employed for lattices of sites. We first 
shall present results on evaporation kinetics, with wave- 
length M ranging from 32 to 512. 

In the DOF phase, K = and, mostly, ksT/L = 0.7, 
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FIG. 4: Time, measured in Monte Carlo steps per site, di- 
vided by M^, dependence of the amplitude Hm of sinusoidal 
gratings for various system sizes, M, and initial amplitudes 
A, setting K — and kBT/L= 0.7. Evaporation kinetics is 
simulated, averaging over several (ranging from A'^ = 10 for 
M = 512 up to = 48000 for M = 32) realizations. The 
dashed lines correspond to an exponential decay of the am- 
plitude Hm- In the inset, the effective relaxation time, r^ff, 
as a function of Hm is shown for M — 64, A = 10. 

the profile of the gratings is quite close to sinusoidal dur- 
ing the equilibration, showing, however, in the beginning 
a tendency towards a broadening near the extrema. The 
time dependence of the decay of the amplitude, 
is illustrated in Fig. 4, studying profiles with fixed ra- 
tio M/A. From this figure and other simulational data, 
we infer that the amplitude falls off at fairly early times 
approximately exponentially, Hm cx: exp(— f/r), r being 
the relaxation time, in accordance with the prediction 
of MuUinsi^ for the rough phase. In that time regime, 
r is proportional to A/^, also as predictedpi^ However, 
as time goes on, in contrast to the behavior in the rough 
phase, a systematic deviation from the exponential decay 
occurs, with a significantly slower equilibration. When 
fitting there the decay to a power-law, the effective expo- 
nent depends both on time and the wavelength M (fixing 
M/A), becoming lower for larger systems (in the largest 
system we studied, M — 512, we find an effective ex- 
ponent of roughly —1.5 to —2). Note that one might 
be still rather far from a characteristic asymptotic tem- 
poral dependence. For comparison, in the flat phase, 
the equilibration is predicted'^^ to follow, at long times, 
Hm oc t~^^'^. The deviation from the exponential form 
sets in at a larger amplitude for gratings with larger 
wavelength and larger initial amplitude A, see Fig. 4. 
Indeed, that amplitude seems to grow like , with 6 
being about 5/4. Thence, we tend to conclude, that the 
theory of MuUins describes merely a transient behavior 
for gratings of flnite wavelength in the DOF phase. 

By inspection of Monte Carlo conflgurations during the 
relaxation, one observes the decay to proceed layer by 
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layer. At a reference moment of time, the top (and bot- 
tom) terrace of the grating may display a flat disordered 
structure, with the sites being predominantly at the in- 
teger heights n and n — I, yielding an amplitude of 
approximately n — 1/2. The terrace is now bordered by 
two meandering steps, and the amplitude of the profile 
(as it follows from averaging over many realizations) may 
decrease fairly slowly, when time goes on. Eventually the 
two steps may touch each other, giving rise to islands of 
flat disordered structures on top of another flat disor- 
dered terrace of average height n — 3/2. The islands then 
shrink and dissolve quite quickly, associated with a rather 
fast decay of the amplitude. Thence, we expect two time 
scales, due to the meandering of the steps and the dis- 
solving of islands, in accordance with previous simula- 
tional results on the equilibration of gratings in the flat 
phasei2SiS&. Indeed, as exemplified in the inset of Fig. 4, 
the two time scales lead to an oscillatory behavior in the 
effective relaxation time, Teff = —{dhiH„i/dt)^^, as a 
function of the amplitude Hm- The pronounced increase 
of Teff at small amplitudes corresponds to the onset of 
the slower, non-exponential decay, discussed above. In 
the rough phase, these oscillations tend to fade away, for 
obvious reasons. 

In the case of surface diffusion, similar features hold. 
We set if = and ksT/L^ 0.7 in the DOF phase 
{ksT/L — 2.5 in the rough phase), with A ranging from 
3 to 6, and with Af ranging from 24 to 48. The equili- 
bration proceeds on much longer time scales than in the 
evaporation case. First, as it is true for evaporation ki- 
netics as well, the discretized initial configuration relaxes 
to a close-to-sinusoidal profile with disordered fiat struc- 



tures at the extrema. Then, still at early times, the de- 
cay is described well by an exponential time dependence, 
with the relaxation time being proportional to M"*, in 
agreement with MuUins' theory for surface diffusion in 
the rough phase.^"* Eventually, as time goes on, in the 
DOF phase a slowing down sets in, like for evaporation- 
recondensation. A crossover from an exponential decay 
to a slower decay at larger times had been found before 
for surface diffusion in the flat phase, with the crossover 
time depending on temperature^ 

Finally, we studied the flattening of one-dimensional 
wire profiles. In the DOF phase, applying evaporation 
kinetics, the equilibration proceeds again layer by layer, 
driven by step annihilation at the disordered flat top ter- 
races, accompanied by two time scales. Analogous fea- 
tures have been found before for the flat phase^'^. The 
excess mass at the surface, due to the wire, changes with 
time, with the average height tending to approach a half- 
integer value. In contrast, in the rough phase, the excess 
mass is nearly conserved. In the rough phase, the am- 
plitude of the wire approaches a decay Hm oc t~^/^, as 
predicted by the MuUins' theory and already conflrmed 
in previous simulations. To determine the asymptotic de- 
cay law in the flals^ii and DOF phases, presumably wires 
with very large initial width and amplitude are needed, 
being beyond the reach and scope of the present study. 

In summary, the preroughening transition is found to 
be non-universal. Dynamic properties in the DOF phase 
resemble partly ones observed in the flat phase, partly 
ones observed in the rough phase. They deserve to be 
studied in even more detail in the future. 
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